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ABSTRACT 

We investigate the orientation of the axes and angular momentum of dark matter halos with respect 
to their neighboring voids using high resolution N-body cosmological simulations. We find that the 
minor axis of halos tends to be aligned along the line joining the halo with the center of the void, 
and that the major axis is perpendicular to this line. However, we find that the angular momentum 
of halos does not have any particular orientation. These results may provide information about the 
mechanisms whereby the large-scale structure of the Universe affects galaxy formation and cast light 
upon the issue of the orientation of galaxy disks with respect to their host halos. 
Subject headings: galaxies: halos — large-scale structure of the universe - methods: statistical — 
dark matter 



1. INTRODUCTION 

The knowledge of the intrinsic alignments of the an- 
gular momentum of dark matter halos with respect to 
their axes and the alignments of these axes with respect 
to the galaxy disks is of great importance in order to 
understand the processes involved in galaxy formation. 
Moreover, the alignment of halo axes, angular momen- 
tum and galaxy spins with respect to the large-scale 
structure may provide further clues to this end; for exam- 
ple, they could be determinant to understand the origin 
of the angular momentum of galaxies (see e.g. Navarro 
et al. 2004). They also may help to understand sev- 
eral observational results like the distribution of galactic 
satellites (e.g. Zentner et al. 2005, Azzaro et al 2005), 
the formation of galactic warps (e.g. Dekel & Shlosman 
1983; Toomre 1983; Bailin 2004) and the tidal streams 
observed for example in the Milky Way (e.g. Ibata et 
al. 2002). Furthermore, they could be a source of con- 
tamination in weak lensing studies (e.g. Heymans et al. 
2006). 

Dark matter halos are in general triaxial being defined 
by their three semi-axes, a > b > c (Dubinski & Carl- 
berg 1991; Warren et al. 1992; Cole & Lacey 1996; Jing 
& Suto 2002; AUgood et al. 2006). Recently, some ef- 
fort have been devoted to measure their intrinsic align- 
ments with respect to the angular momentum. It has 
been found that the angular momentum of dark matter 
halos tends to be parallel to their minor axis (Bailin & 
Steinmetz 2005; Allgood et al 2006). This situation is 
typical for rotationally supported objects. However, as 
dark matter halos are found to be not rotationally sup- 
ported, this alignment could have a non trivial origin 
(Allgood et al 2006). 

On the other hand, the measurement of the alignments 
of galaxy disks with their host dark matter halos is both 
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theoretical and observationally difficult. Using cosmo- 
logical hydrodynamic simulations of disk galaxy forma- 
tion, Bailin et al. (2005) found that the halo minor axis 
within Q.lRyiT. matches the disk axis. Beyond O.lRmr, 
the halo orientation is uncorrelated with the disk axis. 
Moreover, Hoekstra, Yee & Gladders 2004, using galaxy- 
galaxy lensing measurements find that galaxies (disks 
plus ellipticals combined) are aligned with the major axis 
of their dark matter halo. Other works show that the 
angular momentum of disks and dark matter are not 
perfectly aligned (van den Bosch et al. 2002; Sharma 
& Steinmetz 2005). Hence, the final orientation of the 
galactic disk in a cosmological dark matter halo is still 
an unresolved problem. 

One way of gaining insight about these alignments is 
to study the alignments of both, halos and disks, with re- 
spect to the large-scale structure. The basic idea is that 
the angular momentum of halos and galaxies is generated 
by tidal torques operating on the primordial material des- 
tined to form them (Peebles 1969; White 1984). In this 
scenario, the primordial angular momentum of barions 
and dark matter are aligned and the final galactic angu- 
lar momentum is a 'relic' of the primordial one. Accord- 
ing to this, the current misalignment between galactic 
and dark matter halo angular momentum is due to subse- 
quent tidal torques afi^ecting dark matter halos. However, 
this scenario could be more complex because additionally 
to the primordial tidal torques operating on the material 
which eventually becomes the disk, the mass accretion by 
the disk is a clumpy and stochastic process (see Vitvitska 
etal. 2002). Along this line, Navarro et al. (2004) found, 
using high resolution numerical simulations, that galaxy 
disks are rather well aligned with the angular momentum 
of the surrounding matter at the epoch of turnaround. 
Bailin & Steinmetz (2005) found that the minor axes 
of halos of groups and clusters show a strong tendency 
to lie perpendicular to the direction of filaments. How- 
ever, it has been suggested that although filaments cause 
alignments of cluster-size halos, they have not discernible 
influence on smaller halos (Altay et al. 2006). On the 
other hand, it has recently been found using observa- 
tions from the 2dFGRS and the SDSS, that galactic disk 
planes tend to be perpendicular to the walls of cosmic 
voids (Trujillo, Carretero & Patiri 2006). 
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In this letter we study, using high resolution N-body 
simulations, the orientation of axes and angular mo- 
mentum of galactic size dark matter halos embedded in 
the walls of large voids, those with radius larger than 
10/i~^Mpc. We will discuss the observational implica- 
tions of our results. 

2. NUMERICAL SIMULATIONS AND METHODOLOGY 

2.1. N-body simulations 

We have used for our analysis a high resolution N- 
body numerical simulation of 512^ particles in a box of 
120/i~^Mpc on a side assuming the concordance ACDM 
cosmological model (Ho = 0.3, JIa = 0.7 and Hq = 
lOO/ikm/s/Mpc with h = 0.7). The simulation was done 
with the Adaptive Refinement Tree code (ART, Kravtsov 
et al. 1997). The force resolution was 1.8/i~^kpc co- 
moving and the mass of a particle 1.1 x 10^/i~^Mpc. 
Dark matter halos are identified in the simulation by 
the Bound Density Maxima algorithm (BDM, Klypin & 
Holtzman 1997; Klypin et al. 1999). For more details 
about this simulation see Table 1 (I/I2O0.9) in Allgood et 
al. (2006). There is a total of 181000 halos in the whole 
simulation box at z = 0. However, we limited our anal- 
ysis to well-resolved halos with more than 1000 particles 
inside Rmr to be able to analyze accurately their internal 
axes and angular momentum orientations. 

2.2. Searching for voids and halo sample 

In order to study the alignments of dark matter halos 
with their surrounding large-scale structure, we searched 
for large voids in the numerical simulation box using 
the HE void finder (Patiri et al. 2006). The HE Void 
Finder is conceptually simple, as it searches for the max- 
imal spheres with radius larger than a given value that 
are empty of halos with mass larger than some specified 
value. This code is designed to be accurate and computa- 
tionally efficient in finding the largest voids in any point 
distributions like galaxy samples or dark matter halos 
samples. To search for these voids, the code first gen- 
erates a sphere of fixed radius randomly located within 
the sample. We check if this sphere is empty of previous 
selected halos. If so, we increase the size of the sphere 
until its surface reaches the four nearest halos. This is 
our void. For a full description of the algorithm see Patiri 
et al. (2006). 

In the present work, we have searched for voids de- 
fined by dark matter halos with mass larger than 1.0 x 
IO^^/i-^Mq. We choose this mass in order to compare 
with the observational results found by Trujillo, Car- 
retero & Patiri (2006), i.e. we select the halo mass so 
that their number density is equal to the number density 
of galaxies found in the observed galaxy samples (galax- 
ies brighter than Mbj = —19.4 -I- 5logh). With these def- 
initions, we found a total of 60 voids with radius larger 
than 10/i~^Mpc in the whole simulation box. Once we 
have obtained the voids, we selected the halos located 
within shells of 4/i^^Mpc thick beyond the surface of 
the void. Also, as we are interested in the orientations 
of galactic size halos with respect to their neighbouring 
large-scale structure, we restricted the mass of the halos 
to be smaller than 1.0 x IO^^H'^Mq. We detect a total 
of 1729 halos in this mass range that are located in the 
shells of large voids. For each of these halos we obtain 



their principal axes (a, b, c) and calculate their angular 
momentum, in order to measure their orientations with 
respect to the void centers (see next section). 

2.3. Measuring halo orientations 

We follow the method presented in Allgood et al.(2006) 
to fit our halos to ellipsoids. First of all, we calculate the 

reduced inertial tensor for all the particles inside i?vir 
/r^, where the sum is over every parti- 
cle inside a sphere of radius i?vir around the center of 
the halo. The halo center and i?vir are given by our 
BDM halo finder. The coordinates here are Cartesian 
and the denominator (squared distance to halo center) 
makes the eigenvalues of the tensor, Ai > A2 > A3 di- 
mensionless and non-sensitive to particle distances. The 
axis lengths a > b > c are proportional to the square 
root of the eigenvalues. After this first step, we iterate 
this procedure to get more and more approximate ellip- 
soids that fit our dark matter distribution. To achieve 
it, we take the eigenvalues that we have just calculated 
and we normalize them so that a = _Rvir and the ra- 
tios c/a, b/a are preserved. We then take into account 
every particle inside the ellipsoid (the ones that obey 

i^f + (¥) + i^f < 1' witti a, S, c being the imit 

eigenvectors), and calculate the reduced inertial tensor, 
defined now as: 

I,.=J2^ r^ = aVxya^+yyb^ + zyc^ (1) 

Its diagonalization results in the new eigenvalues and 
axis lengths. This method is repeated until convergence 
is achieved in every axis, with a tolerance e = 5 x 10~^. 
We stop the fit after just 7 iterations, as the convergence 
is usually very fast. The halos that could not be fitted 
to ellipsoids (~ 8%) are therefore rejected. 

Once the ellipsoid axes are measured, we calculate in a 
straightforward way the angle between the vector joining 
the center of the halo to the center of the void R = 
RR = Rhaio — -Rvoid, and the principal halo axes with 
a = arccos |.R • a| f3 = arccos \R- b\ 7 = arccos \R ■ c\. 
The la error for the direction of the axes is about 8° (see 
equation 5 in Bailin & Steinmetz 2004). 

As the particles inside the ellipsoid are known after the 
iterative method converges, we can calculate the total 
angular momentum vector of the halo as: 

L = ^mfx V (2) 

where f is again the particle position referred to the 
center of the halo, and v is the particle velocity. Note 
that we are approximating this center with the real center 
of mass of the ellipsoid, otherwise there would be non- 
negligible additional terms to take into account. Then 
the angle 6 between L and R is calculated as follows: 

,R-L, 

6 = arccos (3) 
RL 

There are several sources of noise when calculating 
these angles, so it is not wise to take every halo we have in 
our sample. In particular, we reject the halos whose axes 
are determined to within an uncertainty greater than 
0.2 rad ~ 11°. Moreover, we reject halos with degenerate 
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axes, as the direction of these axes is ill-defined and arbi- 
trary within a plane (see Jing & Suto 2002). Therefore, 
we just take halos with axis ratios lower than 0.95, in 
order to avoid such degeneracy. After all, we have 1187 
halos left for our analysis. 

We also checked if our halos have the expected prop- 
erties reported in the literature for their axis ratios and 
the intrinsic alignment of the angular momentum with 
the direction of principal axes. No discrepancies were 
found. The mean axis ratios in our mass range are 
c/a = 0.64±0.14, h/a = 0.78±0.15, and c/b = 0.82±0.09 
in agreement with Jing & Suto (2002), Bailin & Stein- 
metz (2005), and AUgood et al. (2006). The median 
direction cosines of the angular momentum with the prin- 
cipal axes are L-a = 0.32, L-b = 0.48, and L-c = 0.69, in 
agreement with the results in Bailin & Steinmetz (2005). 
About 75 % of halos have an angle between their angu- 
lar momentum and their major axes greater than 60° as 
shown in Gottlober & Turchaninov (2006). 

3. RESULTS AND CONCLUSIONS 

In Figure 1 we show the orientations of the three halo 
axes and angular momentum with the cosmic voids. In 
the top-left panel we show the probability distribution 
of the a angle which is defined between the major axis 
of the dark matter halos and the vector pointing to the 
center of the void (black dots). Similarly, in the top- 
right and bottom-left panels we present the angles /3 and 
7 which are the ones between the middle and minor axis 
with the direction to the center of the void respectively. 
Finally, in the bottom-right panel we show the angle 9 
between the halo angular momentum and the direction to 
the center of the void. The solid line corresponds to the 
probability distribution of randomly distributed angles, 
which is a sine function for geometrical reasons. 

In these plots we clearly see that halos are not ran- 
domly oriented within the shells of voids. The major 
axis of halos tends to lie along the plane perpendicular 
to the radial direction, i.e. there is an excess of halos 
with a angles near 90 degrees. On the contrary, the mi- 
nor axis tend to point to the direction of the center of the 
void. There is a marginal orientation of the middle axes 
toward the plane perpendicular to the radial direction. 
However, we do not see any particular orientation of the 
halo angular momentum with respect to the large-scale 
structure, confirming recent results obtained also from N- 
Body simulations using similar methodology (Heymans 
et al. 2006). 

To test the robustness of our results we have conducted 
several tests. First of all, we have analyzed the orienta- 
tions of dark matter matter halos selected as described 
above, but in this case we used "fake" voids, i.e. we ran- 
domly placed spheres of fixed radius over the sample and 
calculated the angles between the axis (and angular mo- 
mentum) of the dark matter halos located in the shells of 
these "fake" voids and the direction to their center. As 
expected, we found that the angles of the halos chosen 
to be in the shells of random spheres are uniformly dis- 
tributed, i.e. they follow a sine distribution. Moreover, 
we applied two difi'erent statistical tests: the average of 



the cosine test and a standard Kolmogorov-Smirnov test. 
If the angles are randomly distributed, the average of the 
cosine of these angles is 0.5. So that, the bigger is the 
deviation from this value the larger will be the effect. In 
Table 1 we show the results of the statistical tests. In 
the upper part of the table we show the results for the 
mean cosine (column 1) and the probabilities that the re- 
sults be compatible with isotropy as given by the mean 
cosine and the KS tests (columns 2 and 3) for halos lo- 
cated in the shells of voids and, in the bottom part the 
same results but for halos inside the shells of the random 
spheres (i.e. the 'fake' voids). From these results, we 
conchide that we can reject with high confidence (more 
than 7 sigmas) the null hypothesis (i.e. that the axis of 
dark matter halos are randomly distributed) for the ha- 
los in the shells of voids. As expected, the results of the 
statistical test for the halos located in the shells of the 
random spheres show that the orientation of these halos 
is compatible with the random distribution. 

On the other hand, we used our data to gain some 
insight into the nature of the internal alignment of the 
minor axis with angular momentum found in previous 
works (e.g. Bailin & Steinmetz 2004). We selected those 
halos where both directions form an angle smaller that 
20°. For these halos, we found that the alignment be- 
tween the minor axis and the radial direction of the void 
is weaker than for randomly chosen halos, being even 
compatible with a random distribution. 

Comparing these results with the observed tendency 
for the spin of spiral galaxies to point along the plane 
perpendicular to the direction of the center of the void 
found by Trujillo, Carretero & Patiri (2006), three sce- 
narios are possible: 1) The non- isotropy of the orienta- 
tion of the spin of the galaxies is entirely induced by the 
halo alignments, 2) both alignments are two efi'ects with 
common origin (namely, the tidal field in the shells of 
voids), or 3) a combination of the previous ones. In the 
first scenario the probability distribution for the angle 
between the major axis of the host halo and the spin 
of their galaxy is independent of the angle between the 
spin and the radial direction of the void. In this case, this 
probability distribution may be uniquely determined by 
comparing the distribution of the angle between the ma- 
jor axis and the radial direction with that of the angle 
between the spin of the galaxies and the radial direction. 
This would imply that the spin is very closely aligned 
with the major axis of the host halo. However, recent 
hydrodynamical simulations (Bailin et al. 2005) does 
not show this effect (they found that the spin of galaxies 
is uncorrelated with the halo). Therefore, the other two 
scenarios seem more plausible. 
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TABLE 1 
Results of the statistical tests 



Angle < cos > < cos > Prob. KS Prob. 



Halos in walls of voids 

a 0.459 1.13 x IQ-^ 5.4 x lO"'^ 

/3 0.481 7.58 X IQ-^ 2.0 x 10"^ 

7 0.556 9.19 x lO'^^ 9.4 x 10"" 

e 0.509 0.19 0.12 



Halos in walls of random spheres 

a 0.496 048 023 

13 0.509 8.90 X 10"^ 7.80 x 10"^ 

7 0.492 0.19 0.10 

e 0.502 0.69 0.57 
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Fig. 1. — Probability density distribution for the angles between the axes (and angular momentum) of the dark matter halos and the 
vectors pointing to the center of the voids. Top-left panel: ct angle for the major axis. Top-right panel: /3 angle for the middle axis. 
Bottom-left panel: 7 angle for the minor axis. Bottom-right panel: d angle for the angular momentum. The error bar on each bin is the 
Poissonian error. The solid line represents randomly distributed angles (i.e. a sine distribution) 



